Examination of the Gunion-Bertsch formula for soft gluon radiation 
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The spectrum of emitted gluons from the process g + g— > g + g + g has been evaluated by re- 
laxing some of the approximations used in earlier works. The difference in the results from earlier 
calculations have been pointed out. The formula obtained in the present work has been applied to 
estimate physical quantities like equilibration rate of gluons and the energy loss of fast gluon in the 
gluonic plasma. 
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The radiation of soft gluons from the generic partonic 
(V) processes V1+V2 -+P3+ 7^+gluon play a crucial 
role in the study of quark gluon plasma (QGP) expected 
to be formed in heavy ion collisions (HIC) at ultra- 
relativistic energies. The number nonconserving process, 
g + g— >g + g + g has drawn particular attention in view 
of its importance for the (i) chemical equilibration in 
the deconfined phase of quarks and gluons [H-dj], (ii) 
energy loss of fast gluons propagating through gluonic 
plasma [3, (hi) evaluation of transport coefficients 
of the gluonic plasma ||6^9j, etc. The Gunion-Bertsch 
(GB) [lOj formula for the spectrum of the radiated glu- 
ons from the processes g + g— >g + g + g has been widely 
used in many of these calculations. Recently, attempts 
have been made to revisit/generalize the GB formula 
[ill [I2I ]. The main purpose of the present work is to 
hnd the corrections to the GB formula by relaxing some 
of the approximations previously adopted. Then, we in- 
spect the effects of the correction terms relative to the 
GB formula on some physical quantities like the equili- 
bration time and the energy loss of fast gluons in a glu- 
onic plasma. The results obtained here can be extended 
to other partonic processes such asg + g— s-g + g + g, 
q + g — » q + g + g, etc., where q stands for quark. 

We consider the process g(/ci)+g(fc 2 ) — > g(fc 3 )+g(fc 4 ) + 
g(ks). The square of the invariant amplitude for this 
reaction can be written elegantly as [la ]: 
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where 



Af = (fci.fc 2 ) 4 + (fci.fc 3 ) 4 + (fci.fc 4 ) 4 

+ (fcx.fcs) 4 + (fc 2 .fc 3 ) 4 + (fc 2 .fc 4 ) 4 + (fc 2 .fc 5 ) 4 

+ (fc 3 .fc 4 ) 4 + (h.k 5 ) 4 + (fc 4 .fc 5 ) 4 ; (2) 
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and 



V = (fc 1 .fc 2 )(fc 1 .fc 3 )(fc 1 .fc 4 )(fc 1 .fc 5 )(fc 2 .fc 3 ) 
x (fc 2 .fc 4 )(fc 2 .fc 5 )(/c 3 .fc 4 )(fc 3 .fc 5 )(fc 4 .fc 5 ), 

(ijklm) = (ki.kj)(kj.kk)(kk-ki)(ki.k m )(k m .ki) 
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N c (= 3) is the number of colors, g = -y/47ra s is the color 
charge, and a s is the strong coupling. 

The quantity, |M gg _>. ggg | 2 after simplification can be 
written as (see Appendix): 
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where |M gg ^ gg | 2 

{k\ — fc 3 ) 2 , u = (k\ — fc 4 ) 2 , k± is the transverse mo- 
mentum of the radiated gluon. The 0(kj 1 ) and 0{k Q L ) 
terms appearing in Eq.[5]were absent in Ref. Hence- 
forth these two terms will be called the correction terms. 
We will demonstrate that the contributions from these 
terms are non-negligible and will have crucial importance 
for the phenomenology of heavy ion collisions at ultra- 
relativistic energies. 

While the details for the derivation of the Eq. [5] is 
given in the Appendix, we would like to check the effects 
of the correction terms in Eq. [S] to physical quantities 
like equilibration time of gluons and energy loss of fast 
gluons propagating through a gluonic fluid. 

Let us first discuss the role of the correction terms 
in the equilibrium time scale of gluons. The processes 
g + gf^g + g, g + g«g + g + g, g + gog + g + g + g 
etc are responsible for maintaining equilibration in the 
system. In the present work we would like to estimate 
the contribution from the correction terms of the process 
g + g-H-g + g + gas mentioned above. We evaluate the 
equilibration rate for this process (with s = 18T 2 as in 
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with and without the correction terms. The equili- 
bration rates obtained from the spectra of Refs. [ll|, 12 1 



2 



and present work normalized by the GB spectra (r^) are 
displayed in Fig. [TJ We observe that the equilibration 
rate obtained with the correction terms is smaller com- 
pared to the scenario when the corrections are neglected 
for the entire range of temperature under consideration. 
This indicates that the contribution from the correction 
terms will enhance the equilibration in the gluonic sys- 
tem. 

Before going into estimating the energy loss of fast 
partons moving through a gluonic fluid, some clarification 
are in order here. The energy loss of fast partons in QGP 
has been evaluated rigorously (see [H], [H[ f° r review) . In 
the present work the aim is not to achieve the same level 
of rigor but to check the effects of the correction terms 
discussed before on the energy loss of fast partons. 

The energy loss of high-energy partons propagating 
through QGP is a field of high contemporary interest. 
Experimentally the energy dissipation has been measured 
through the suppression of the transverse momentum 
(Pt) distribution of hadrons produced in Au+Au col- 
lision relative to the binary scaled p+p interaction at 
the same center of mass energy. The nature of the sup- 
pression may be used as a tool for diagnosis of QGP 
formation in nuclear collisions at Relativistic Heavy Ion 
Collider (RHIC) and Large Hadron Collider (LHC). The 
two most common mechanisms for the energy loss are 
elastic and inelastic or radiative processes. Among vari- 
ous in-elastic processes involving quarks and gluons, the 
process: g + g— >g + g + g plays a major role. 

In the hard processes the gluons produced with virtu- 
ality ~ Q 2 is highly off-shell because during the collision 
process the color field of the gluon is stripped off. There- 
fore, the highly virtual gluon will develop a dead cone 
- which will results in radiative suppression. Following 
the procedure outlined in Ref [l6| we evaluate the en- 
ergy loss of gluon, AE(L) as a function of length, L. 
The results obtained with GB and present gluon spectra 
are displayed in Fig. [2] We observe that the difference in 
AE(L) between the GB spectra and the spectra derived 
in the present work (with corrections described above) is 
small at lower L but grows up with increase in L. 

The spectrum of the radiated gluon in medium, derived 
by using the ratio of the amplitude square of the radiative 
process gg — > ggg to that of the elastic process, gg — » gg 
is given by, 
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where r\ is the rapidity of the radiated gluon, the sub- 
script GB has been used to indicate the gluon s pec trum 
obtained using the approximation considered in [10J (see 



also [13]) which is given by, 
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where tud 



2Tta s (T){C A + ^f)T/3, is the thermal 



mass of the gluon [18|, Np is the number of flavors con- 
tributing in the gluon self-energy loop, Ca — 3 is the 
Casimir invariant for the SU(3) adjoint representation, 
a s is the temperature-dependent strong coupling [l9| and 
q± is the transverse momentum transfer. The thermal 
mass in the denominator of Eq.[7]has been introduced to 
shield the infrared divergence arising from the massless 
intermediary gluon exchange. 

The energy loss of a gluon passing through a gluonic 
medium can now be calculated using the gluon spectrum 
of Eq. [6] The energy loss per collision can now be esti- 
mated as: 
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where fco = k±coshr) is the energy of the radiated gluon 
and tf is the formation time of the gluon. The first 9- 
function in Eq. |8l involving A -1 or interaction time, is 
introduced for the Landau-Pomeranchuk-Migdal (LPM) 
effect. The LPM effect imposes restriction on the phase- 
space of the radiated gluon, it must have rp(= coshr]/k±) 
less than the mean free time, A -1 . The gluon can be 
emitted for time scales larger than tf [![. The second 
^-function sets the upper limit for the energy of the ra- 
diated gluon. 

To proceed further, we replace q\ by its average value 
evaluated as follows: 
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Ci is 9/4, 1 and 4/9 for gg, qg, and qq scattering. (q±) 
is then obtained as, 
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For y/s —> oo, i.e. in the hi gh- e nergy limit one can make 
the replacement t ~ — q\ [20|. In contrast to previous 
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works [ll|, the value of s is taken as s ~ 6£T [4] 
allowing the possibility for the incident gluon to remain 
out of thermal equilibrium. With all the above ingredi- 
ents we are now ready to evaluate energy loss (dE/dx) 
of a fast gluon in a gluonic medium as follows: 
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The interaction rate, A has been evaluated by using the 
procedure similar to [2lj . We need to implement now 
the radiative suppression of gluons due to its possible 
off-shellness (Fig. [2]). The energy loss (AE) of off-shell 
gluons of energy 15 GeV moving through a gluonic fluid 
of dimension 4 fm for the gluon spectrum derived in 
Refs. [llj], [HI and in the present work [obtained from 
Eqs. [fO [8] and [13] normalized to that resulting from 
GB approximations (AEr) are displayed in Fig. [31 We 
observe that with the correction terms the value of AE 
is enhanced by about 40% and 20% for T 300 MeV and 
400 MeV, respectively, compared to the AE obtained 
from the spectra of Refs. [ill ] and [TJ|]. Such differences 
may have important consequences on the heavy-ion phe- 
nomenology at RHIC and LHC collision energies. 

In summary we have evaluated the spectrum of the 
emitted gluon from the processes g + g— > g + g + g by 
relaxing some of the approximations used in recent calcu- 
lations. The results derived in the present calculation has 
been applied to evaluate some physical quantities, e.g. 
energy loss of energetic partons in HIC and time scale for 
equilibration in gluonic plasma. Results obtained in the 
present work have been compared with the earlier calcu- 
lations [ll|, [TH . We find that the contributions from the 
previously neglected terms 0(kj 1 ) and O(k L ) in the ma- 
trix element play crucial role in heavy-ion phenomenol- 
ogy. The gluon spectrum obtained in Eq. [6] can be used 
for other 2 — > 2 + g partonic processes also. 
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FIG. 1: (Color online) Temperature variation of the ratio 
of the equilibration rate (inverse of the time scale) obtained 
in the present work (solid line), Ref. [Til ] (dashed line), and 
P2] (dot-dashed) normalized by the GB value for the process 




FIG. 2: Energy loss of a 15 GeV gluon in vacuum as a function 
of path length. 
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FIG. 3: (Color online) Temperature variation of AEr of a 15 
GeV gluon moving through a gluonic heat bath of dimension 4 
fm. Solid (dashed) line indicates result for the gluon spectrum 
obtained in the present work ( [Til]). The dot-dashed line 
stands for the results for the gluon spectrum of [L2|. The 
AE normalized by the corresponding value obtained from GB 
approximation and the resulting quantity is called AEr. 



Appendix 

In this appendix we derive Eq. [3] for the square of the 
invariant amplitude for the radiative process, gg — » ggg 
upto orders O(kj_ ) and 0(^). Consider the reaction 



;(fci) + g(*a) -> g(*s) + g{ki) + g(h 



(14) 



where fcs is the four-momentum of the radiated gluon. 
The Mandelstam variables for the above process are de- 
fined as 



s = (h + fc 2 ) 2 , t = (k 1 - k 3 ) 2 
u = (ki — fc 4 ) 2 , s' = (k 3 + fc 4 ) 2 



t' = (fc 2 -fc 4 ) 2 
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Because gluons massless we can write 
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We also have the relations 

s + t + u , , s + t'+u' 

ki-ks = ^ ' fe2 - fe5 = 2 

s + t' + u , , s + t + u' 

k 3 -k 5 = , fc 4 .fc 5 = 7. • 



(16) 



(17) 



u(s + u> + t')(s +u + t') + ss'' 

l - [J_ 

t(s + u' + t'){s +u' + 1) uu' 

1 , 



and Af can now be written as 

TV = ^ + t i + u 4 + s' 4 + t' 



/4 



(23) 



+ {s + t + uf + (s + t' + u') A + {s + t' + uf 
+ {s + t + u') 4 }. (24) 



For soft gluon emission, 

s + t + u + s' + t' + u' = 0. 



(18) 



The matrix element square of the radiative process 
gg -> ggg is given by jlj| 



|M g 



2 9 N 2 ~IV 
x [(12345) + (12354) + (12435) 
+ (12453) + (12534) + (12543) 
+ (13245) + (13254) + (13425) 
+ (13524) + (14235) + (14325)], (19) 



where N c is the number of colors, g = V '47ra s is the 
strong coupling, 

Af = (fc 1 .fc 2 ) 4 + (fc 1 .fc 3 ) 4 + (fc 1 .fc 4 ) 4 

+ (fcx.fcg) 4 + (fc 2 .fc 3 ) 4 + (fc 2 .fc 4 ) 4 + (fc 2 .fc 5 ) 4 

+ (k 3 .h) 4 + (k 3 .k 5 ) 4 + (fc 4 .fc 5 ) 4 , (20) 

V = (fa.fc)(ki.fe)(ki.*a)(*i.*B)(fe.*&) 

x (k2.k4){k 2 .k 5 ){k 3 .k i )(k 3 .k 5 )(k i .k 5 ), (21) 

and 

(ijklm) = {k i .k J ){k J .k k ){k k .ki){ki.k m ){k m .k i ). (22) 
Simplifying Eq. [THlwe get, 



For a soft gluon emission (k$ — > 0) s —> s', t — >• t' , 
u — > u' . We can express the transverse momentum of 
the emitted gluon as 

k 2 ± = 4(fc 1 .fc 5 )(fc 2 .fc 5 )/s 

= (s + t + u)(s + t' + u') 
= (s + t + u) 2 /s 

(25) 

Using Eqs. [23l [24] and [25l the square of the matrix ele- 
ment can be written as 
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where the subscript GB stands for the approximation 
used by Gunion and Bertsch [l0|. For the elastic process, 



l^gg->gg| 2 = V 
On simplifying Eq. [26] we obtain, 
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In the proposed kinematic limit we set terms which are 
linear in k± to zero and keep terms C(fc° ), ©(fcj 1 ), and 

C(fcJ 2 ) in |A^| 2 ^ ggg - We also neglect terms O(^) and 
higher order in the matrix element. To proceed further 
one needs to express the Mandelstam variable u in terms 
of s, t, and k± by using the following relation: 
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The binomial expansion of [1 — — -)] 1 converges if 
f*i _ *) < i, F or the kinematic limit mentioned above 

i.e. for k± —> and keeping terms upto 0(^), the in- 
equality — |) < 1 is satisfied. We have checked that 



terms beyond — -|) 5 in the expression of i are not re- 
quired for the kinematic limit under consideration. With 
all these we get, 
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For the assumed kinematic conditions u 4 /s can be ex- 
pressed as follows: 
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Similarly, 
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Putting Eqs. |3fl]to [Him [28] we get, 
l M ! 2 gg-ggg = 12 3 2 |Af g «| 2 
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The terms ©(fcj 1 ) and O(k x ) contribute to the energy 
loss of the gluons in a gluonic plasma and hence are 
important for heavy-ion phenomenology at RHIC and 
LHC energies. These terms were absent in the previous 
work [I2| (also in [HI) 
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